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ABSTRACT 

The enigmatic outflows of the Orion-KL region have raised discussions about their potential driving 
sources for several decades. Here, we present C 18 0(2-l) observations combined from the Submillimcter 
Array and the IRAM30m telescope. The molecular gas is associated on large scales with the famous 
northwest-southeast high-velocity outflow whereas the high-velocity gas on small spatial scales traces 
back to the recently identified submm source SMAl. Therefore, we infer that SMA1 may host the 
driving source of this outflow. Based on the previously published thermal and maser SiO data, source 
/ is the prime candidate to drive the northeast-southwest low-velocity outflow. The source SMAl is 
peculiar because it is only detected in several submm wavelength bands but neither in the infrared nor 
cm regime. We discuss that it may be a very young intermediate- to high-mass protostar. The estimated 
outflow masses are high whereas the dynamical time-scale of the outflow is short of the order 10 3 yrs. 

Subject headings: stars: formation 
outflows 

1. INTRODUCTION 



stars: early-type - stars: individual (Orion-KL) - ISM: jets and 



The explosive large-scale northwest-southeast outflow 
in the Orion-KL region is one of the best studied but 
probably also least understood outflows in the sky. No 
other region shows that many H2 bow shocks distributed 
in a fan-like fashion with red- and blue-shifted emis- 
sion features on both sides of the c entral region (e.g., 
ISchultz et aLlll999t iNissen et al.ll2007ft . Furthermore, the 
molecular emission as traced in CO shows extremely broad 
line- wings up to > 50kms _1 , but the morphological 
structure has so far been uncle ar and difficult to match 
closely to the H 2 outflow (e. g., IChernin k Wright) fl996: 
Rodrigu ez-Franco et a l. 1999). In addition to this out- 
flow, also known as the high- velocity outflow, there exists a 
second outflow in the northeast-southwest direction which 
was first detected by H2O maser emission an d which is of- 
ten labeled as the low-velocit y outflow (e.g.. iGenzel et all 
119811: iGenzel fc Stutzkilfl989ft . 

The major question is which sources drive the differ- 
ent outflows. While originally the infrared source IRc2 
was a suspect, it was soon realized that one of the prime 
powering sources of the region is the very close-by radio 
source J riMenten k Reidl 1 1995ft . This source is the site of 
vibrationally excited SiO maser emission, and the orien- 
tation of these masers allowed different interpretations of 
the associated outflow direction. While early work favored 
source I t o drive the large-scale northwest-southea st out- 
flow (e.g- IWright et alJll995t fGreenhill et al.lll998ft . more 
recent work on the maser data indicates that it rather 
belongs to the low-v elocity northeast-southwest outflow 
(|Greenhill et al.l l2003) . The latter interpretation finds ad- 
ditional support in thermal SiO observations that show an 
elongated structure centered on sourc e I with an orien- 
tation also in the n ortheast-southwest (jBlake et al.l Fl996: 
iBeuther et all 120051 see also Fig. [TJ). Since thermal SiO 
emission is usually attributed to shocks within outflows 

1 Max Planck Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg. 

2 Department of Physics and Astronomy, University of Aarhus, 8000 Aarhus C, Denmark 

1 



(e.g.. ISchilke et aT]|l997t ). the combination of the thermal 
and maser SiO emission strongly favors source / as the 
driving source of the low- velocity northeast-southwest out- 
flow. Source n has also been discussed as a possible driving 
source of the northeast-southwest low- velocity outflow. It 
shows an elongation along position angle 1 00-130° in the 
mid-infrared which could be due to a disk (jShuping et al.l 
[2001 iGreenhill et all 12004) . and a roughly perpendicu- 
lar elongation in the radio at 8.4 GHz wh ich may be 
caused by an outflow (|Menten fc Reidl f"l995ft . The rela- 
tionship between sources I and n in driving the northeast- 
southwest out flow(s ) is unknown (for a discussion see, e.g., 
IShuping et aHl2004t INissen et alJl2007ft . 

That leaves as an open question which source is the 
culprit for driving t he large-scale high- velocity northwest- 
southeast outflow? Ballv (2008j) recently suggested that 



this outflow may be powered by the dynamical decay of 
a massive star system, but other scena rios are possible 
as we ll. Based on low-resolution data, Ide Vicente et al.l 
(2002) proposed the existence of another luminous source 
approximately 2" south of source I. This source was un- 
ambiguously identified in the first sub-arcsecond resolution 
865 /im submm continuum image obtained with the Sub- 
millimeter Array (S MA), and henceforth labeled SMAl 
(|Beuther et al.ll2004 see also Fig. [J). Follow-up SMA ob- 
servations at 440 ^ m detected the sourc e at these shorter 
wavelengths as well (|Beuther et al.ll2~006ft . In this work we 
show that the high- velocity g traced by C 18 0, appar- 
ently traces back to this relatively unknown source rather 
than to sources I or n. 

2. OBSERVATIONS 

Orion-KL was observed with the SMA in the compact 
configuration simultaneously at 690 GHz and 230 GHz 
on 2005 February 19th. To accomplish the high- 
frequency observations the weather conditions were excel- 
lent with zenith opacities r(230GHz) between 0.03 and 
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0.04 throughout the night. The phase center was the 
nominal position of source I R.A.(J2000) 05h35ml4.50s 
and decl.(J2000) 05°22'30."45, and the w lsr was ~5kms" 1 . 
The double-sideband receivers covered the frequency 
ranges 218.85 to 220.85 GHz and 228.85 to 230.85 GHz. 
The phase was calibrated with regularly interleaved ob- 
servations of the quasar 0607-157, flux and bandpass cal- 
ibration were performed with Callisto observations. The 
high-frequency d ata along with more o bservational details 
are published in iBeuther et al.l (|2006f ) , here we are only 
interested in a small subset of the low-frequency data. 
Complimentary C 18 0(2-1) data were obtained with the 9 
pixel HERA array on the IRAM 30 m telescope in Septem- 
ber 2007. The observations covered a region of 1' x 1' 
around source / and were conducted in the on-the-fly mode 
with on-source integration times per position of ~3.7secs. 
The OFF-source reference position was 15' east of source 
/. The single-dish data were converted to visibilities and 
subsequently processed with the SMA data within the 
MIRIAD package using standard tools (e.g., UVMODEL). 
The weighting of the two datasets was chosen to recover 
large-scale emission but maintain at the same time the 
high spatial resolution to resolve small-scale structure. 
The synthesized beam of the combined data is 4.2" x 3.3", 
and the la rms per 5kms _1 channel is ~50 mJy beam -1 . 

3. RESULTS 

The whole low-frequency SMA dataset will be presented 
in a forthcoming publication by Nissen et al. (in prep.), 
here we are only interested in the C 18 0(2-l) emission. Al- 
though the SMA and IRAM 30 m observations also covered 
the more abundant isotopologues 12 CO(2-l) and 13 CO(2- 
1), in Orion-KL it turns out that because of their high op- 
tical depth these two lines strongly suffer from confusion 
problems with the ambient cloud complicating any inter- 
pretation of the data. In contrast to that, the rarer iso- 
topologue C 18 0(2-l) penetrates the ambient cloud more 
deeply and hence traces the outflow associated gas far bet- 
ter, allowing us a deeper analysis of the outflow structure, 
its mass flow rate, kinematics and energetics, as well as 
pinpointing the likely driving source. Fig. [5] presents the 
combined SMA+30m C 18 0(2-l) channel map. Exclud- 
ing the channels with strongest emission closest to the vi sr 
(approximately four channels between and 15kms _1 ), 
two main features can be distinguished in the remaining 
channels: 

(1) The large-scale gas emission exhibits a V-shaped 
cone-like open structure to the northwest, particularly 
prominent in the red-shifted channels > 20kms~ 1 but 
also visible in the blue-shifted channels between —20 and 
— lOkms -1 . This cone- like structure approximately en- 
compasses the main features of shocked strong H2 emis- 
sion in this direction (Fig. [3]). Finding both blue- and 
red-shifted gas to the northwest indicates a wide opening 
angle of the flow where the main outflow axis is not too 
far from the plane of the sky. The less pronounced emis- 
sion to the southea s t was already previously observed, and 
[Ch ernin fc Wrightl (|1996l ) argue that the outflow may be 
blocked in this direction by the hot core0. It is interesting 
to note that the main southeastern emission features are 

3 The hot core traced in NH3 {Wilson et al. 2000) actually extends fur 
in Fig. [T] The proposed blocking of the outflow by the hot core may 
emission. 



well correlated with a lack of H2 emission indicating par- 
ticularly high column densities in this direction. Further- 
more, we find high-veloc ity gas toward the southwest and 
northeast (Fig. [3]). Since iNissen et all <f2007h find H 2 emis- 
sion with velocities up to 40kms~ 1 toward the southwest 
as well, it is likely that this high- velocity gas is associated 
with the low-velocity outflow. 

(2) We find that the high-velocity gas on the blue- 
and red-shifted side of the spectrum traces back to the 
source SMAl. Although the synthesized beam is larger 
than the spatial separation of source I and SMAl (~ 2"or 
~830AU), the peak positions can be determined with a 
much higher spatial accuracy depending on the signal- 
to-noise ratio S/N. For point sources, one ca n in princi- 
ple a chieve an accuracy of 0.45 HPBW/(S/N) (jReid et alj 
Il988f ). however, this may be less applicable for extended 
emission like C 18 in this case. The channel with the low- 
est S/N ratio still detects the main peak position toward 
SMAl at higher than 15cr values (at -30kms _1 ). For 
a point source, this would imply that the peak positions 
could be determined to better than 0.125". Although the 
C 18 emission is more extended, the shape of the high- 
velocity gas close to SMAl can to first order be described 
by a two-dimensional Gaussian. Therefore, being more 
conservative for the C 18 high-velocity gas, it should still 
be possible to infer the emission peak positions to an ac- 
curacy of better than 1". Based on this, the blue- and 
red-shifted high- velocity gas is spatially mainly associated 
with SMAl, and not with the previously often discussed 
sources / or n. 

The blue- and red-shifted maps shown in Figure [3] (in- 
tegrated from —35 to —20 and from 22.5 to 37.5 kms -1 , 
respectively) can be used to calculate the masses, kine- 
matic and energetic parameters of the outflow following 
Bcu ther et alj ((2002b) (time-scale t, outflow rate M out , 
momentum rate p, and energy E). Assuming optically 
thin C 18 0(2-l) emission at 30 K with a C 18 Q /H 2 conver- 
sion factor of 1.7 x 10~ 7 (|Frerking et al.lli982l ). the derived 
values are given in Table [TJ While the masses have uncer- 
tainties of approximately a factor 5, the kinematic and 
energ etic values are uncertain within an order of magni- 
tude (|Cabrit fc Bertoutlll990l ). The derived masses, out- 
flow and momentum rates are at the upper end of typically 
found massive outflows, whereas the time-sc ale is relatively 
short ( compared with the outflow sample in lBeuther et al] 
I2002bft . 

4. DISCUSSION AND CONCLUSIONS 

The large-scale cone-like morphology of the C 18 out- 
flow approximately following the H 2 bow-shock emission 
in the northwest direction suggests that C 18 and H 2 are 
tracing the same outflow. Furthermore, the spatial associ- 
ation of the C 18 high-velocity emission with the submm 
source SMAl indicates that SMAl may harbor the driv- 
ing source of the enigmatic large-scale Orion-KL outflow. 
Only the combination of the single-dish 30 m data, trac- 
ing the large-scale structure, with the SMA high-spatial- 
resolution observations allows us to see all relevant struc- 
tures and to trace the outflow back to pinpoint its origin. 

her southeast of SMAl than indicated by the hot core continuum peak 
also be responsible for some of the large-scale asymmetries of the H2 
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In combination with the thermal and maser SiO data, 
we suggest that the two driving sources of the two known 
outflows can be identified: (1) SMA1 may host the driving 
source of the large-scale northwest-southeast high- velocity 
outflow, whereas (2) source / is likely the origin of the 
northeast-southwest low- velocity outflow. 

The nature of source I has been subj ect to many in- 
vestigations for more than a decade (e.g.. iMenten &: Reidl 
Il995t IBeuther et all l2006UReid et atll2007f) . however, we 
do not know much about SMA1. So far SMA1 has only 
been detected at submm wavelengths (865 and 440 /mi, 
IBeuther et"aLll2004l I2006D but neit her at shorter infrared 
nor at longer cm wav elengths (e.g.. IMenten fc Reidlll995l : 
iGreenhill et al.l [20041 ) . The submm-only detection is in- 
dicative of a very young age of the source which is sup- 
ported by the short dynamical time-scale of the molec- 
ular outflow. Although deriving a luminosity with just 
two data-points is impossible, based on the large outflow 
masses, it is likely that SMAl is in the process of form- 
ing a massive star. The non-detection at cm wavelengths 
could either be due to quenching of an underlying hyper- 
compact Hn region by very high accretion ra tes and/or 
gravitational tra pping of the ionized gas (e.g., Walmslcv 
□Ml Keto 200!), or the source is still in such a young evo- 
lutionary phase that the central object has not yet enough 
ionizing flux to produce a detectable hypercompact Hn 
region. In the frame work of the recently proposed evolu- 
tionary sequence by IBeuther et alj (|2007l ). it may be an 
intermediate-mass protostar destined to become massive 
in the future. 

Assuming optically thin dust emiss ion at 100 K with 
a dust opacity inde x (3 = 2 (e.g., iHildebrandl [l983; 
IBeuther et al. 2002a!), the spatially resolved submm fluxes 
(|Beuther et al.l 12004120061) correspond to a gas mass anc. 



column density of ~0.1M e and 3 x 10 24 , respectively. 
While the large column densities confirm that this source 
is undetectable at near-infrared wavelengths with current 
instrumentation, at first sight the mass app ears small. 
However, as discussed in (|Beuther et al.|[2004l ). more than 
90% of the total flux is filtered out by the SMA observa- 
tions, and hence the mass estimates are correspondingly 
unreliable. The total luminosity o f the region of about 
10 5 L© (e.g., IMenten k, Reidl fl995f) cannot be produced 
by the sources previo usly discussed in the literature (e.g., 
IGreenhill et al.ll2004D . and it is well possible that SMAl 
contributes a significant fraction of the overall Orion-KL 
luminosity. 

We note that the short dynamical outflow time-scale 
of the order 10 3 yrs is consistent with both recently dis- 
cussed scenarios for the proper motions in this region: 
( 2005) proposed that the recent passage of the BN-object 
has triggered the outfl ow b urst in this region, whereas 
iRodriguez etTal] (|2005h and iGomez et all £2005) suggest 
that a dynamically decaying multiple system is at the root 
of the observed proper motions. Both events are suggested 
to have taken place about 500 yrs ago. Our data do not 
allow to differentiate between these scenarios, however, it 
is interesting that all recent observations are indicative of 
strong activity within the last 10 3 yrs. 
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funded by the Danish Basic Research Foundation and the 
financial support from the Instrument Centre for Danish 
Astrophysics (IDA), funded by the Danish National Sci- 
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Table 1 
Outflow parameter 



M bluo [M Q J 


33 


M rcd [Mq] 


81 


Mtotal [Mq] 


124 


t [yr] 


971 


M out [Moyr- 1 ] 


1.2 x ict 1 


p [Mgkms" 1 ] 


3980 


E [erg] 


1.4 x 10 48 
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Fig. 1. — SiO(8— 7) integrated emission as contour overlay on the grey-scale 865 fim submm continuum e mission obtained with the SMA. This 
figure is an adaption of the 865 /im continuum and SiO(8-7) line data from Beuther et al. (2004, 2005). The grey-scale continuum contours 
start at the 3a level of 0.105 mjy beam -1 and continues in 2<r steps. The SiO emission is contoured from ±10 to ±90% (step ±10%) of the 
peak emission of 2.1Jybeam — 1 with full and dashed contours as positive and negative features (due to missing short spacings), respectively. 
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Fig. 2.— C ls O(2-l) channel map with 5 km s 1 of the combined SMA±30m data. The contours start at the 3<r levels of 150 mjy beam 1 
and continue in 6cr steps for all panels except of the 0, 5 and lOkms -1 panels that follow 12, 12 and 14cr steps for clarity reasons. Negative 
features are draw as dashed contours. The stars mark the positions of the five main continuum source with labels in the top-left panel. The 
synthesized beam is shown at the bottom-left of each panel. 
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Fig. 3. — Overlay of the H2 emission in grey-scale (Nissen et alj|2007h with the blue- and red-shifted C ls O(2-l) emission in color-contours. 
The integration ranges are [—35.0/— 20.0] and [22.5,37.5] kms -1 , respectively. The contours are in 3<r steps of 150 mjy beam -1 . The stars 
mark the positions of the same continuum sources as in Figure [2] the synthesized beam is shown at the top-right of each panel. The large 
circle in the left panel shows the size of the SMA primary beam, whereas the squared box marks the size of the inlay-box shown in the right 
panel. The outflow directions of the high- and low-velocity outflow are sketched in the left panel. 
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